Integral equation of quantum 
stochastic process 
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Abstract 

To describe stochastic quantum processes I propose an integral 
equation of Volterra type which is not generally transformable to any 
differential one. The process is a composition of ordinary quantum 
evolution which admits presence of a quantum bath and reductions 
to pure states. It is proved that generically solutions stabilize asymp- 
totically for t +00 to a universal limit - the projection onto the 
state with maximal available entropy. A number of typical methods 
of finding solutions of the equation are proposed. 
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I Introduction 



In classical evolution of a system-thermal bath it is observed existing of time 
arrow. In quantum case there are a number of attempts to modify equations 
to gain non-invertability of the dynamics. A basis for a formal stochastic 
disturbance is usually a differential equation as the generalized Schrodinger 
equation^ or Lindblat equation, see for example ^. The first modification 
leading to Hugstone equation is a model of single measurement process. Nev- 
ertheless, it realizes one reduction event as a long-time limit and the method 
is not applicable here, where I assume that reductions moments are points 
on the time axis. In the second approach the basis Lindblat equation is only 
a constant-coefficients approximation of complete quantum evolution valid 
for relatively short time regime, compare ^. In ^ it is described a quantum 
stochastic process which avoids formalizing to differential equations at all. 
It consists from mixed quantum evolutions and reductions, where it is as- 
sumed that the reduction skips arise in time moments treated as a Poisson 
stochastic process with a characteristic parameter v. 

In section ^ I propose the new equation governing the evolution of this 
type, which for = is equivalent to the mere quantum equation (^). In 
section ^ the asymptotic theorem is proved. Solutions of some special ver- 



sions of the equation are investigated in section III. In the appendix fulfilling 



of some genericity conditions for the integral equations based on hamiltonian 
evolution are elaborated. 



II Integral equation 

In the article I restrict myself to finite dimensional case. Let a quantum 
evolution of a system in a system-bath pair is defined by 

A„ : M ^ Bin) (1) 
E^-^t = E^l^" = l (2) 

a a 

where 7i = C", B{7i) is the linear, bounded operators set on 7i, a = 
1, . . . , A^, compare also with ^ ^. Now, let = to < • • • < i^j < U+i < ■ ■ ■ 
be a realization of a Poison process with the mean number of events in unit 
intervals equals to z/ > 0. Between two neighboring points in which quantum 
reductions take place tj, tj+i, the evolution is defined by the doubly stochastic 
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matrix (2E) 

M,,it):=J2TrP,A^{t)P,AUt) (3) 

a 

M,,(0) = 1, (4) 

where t := ti+i — ti, an ortho normal basis \i >G T-C,i = l,...,n is chosen and 
Pi := \i >< i\. As one parameter unitary group evolution is determined by 
a self-adjoint operator in quantum bath presence case it is determined by a 
set of operators Bab = bI^. Aab{t), a = ab, are solutions of 

^ab{t) = E BacAcbit) = E Aac{t)Ba (5) 

c c 

with the initial condition Aab{0) = -^^1, where n2 is the dimension of the 
bath Hilbert space, a,b,c = 1,2,..., 

Let M{T) be the average doubly matrix obtaining by summing up all 
realizations of the Poisson process. The result is: 

M{T) = E / dtn dtn^i ... dtiMiT - t„) o . . . 

oM(t2 - h) o M(ti)z/"e-'^^. (6) 

Following a structural similitude to Weinberg- Van Winter equation^ one ob- 
tains from (!]) the integral equation: 

M(t) = e-'^iMiT) + u r M{T - t)M{tY'dt). (7) 

JO 

In the way I have built a dynamical system entirely defined by Bah ^ Bi^H) 
and z/ > by equations (||, |^. The functions A^if) given by are analyti- 
cal, so M(t) as well. Then performing h differentiations of I gain: 

jVf(T)« = e-"^[M('=)(T) + Lfc(T) + z/ C M{T - t)M(t)(^)e"*rft] (8) 

JO 

where are defined by the recurrence formula: 

Lfc+i(T) = z/M(T)M(^)(0) + (1 - z/)Lfc(T) - z/M«(T). (9) 

Generally, the solution of (|^) is always uniquely defined by M(t), an inte- 
grable map from R to the set of general stochastic matrices and given by the 
convergent series (|^). 
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Remark 1. It is appeared that if one takes a different distribution of reduc- 
tions events and still wants to keep an analog of equation it needs to be 
again a Poissonian distribution. Nevertheless, the integral equation general- 
izes itself to the following one: 

M(T) = a{T)M{T) + M{T - t)M{t)b{t, T)dt, (10) 

J 

where a, h are nonnegative, continuous functions such that 

a(T)+ r a{h)b{h,T)dti+ f\{t2)h{t2,ti)h{tuT)dt2dti + . . . = 1 (11) 
Jo JO Jo 

or equivalently 

/ b{t,T)dt = -a{T) + l. (12) 
Jo 

At the moment still if M{t) are stochastic matrices (or 2S) then M{t) are of 
the same type. 



Ill Solving the equation 

The following property of doubly stochastic matrices will be of interest. 

Definition 1. For a doubly stochastic matrix M : M" R" its compression 
c(M) is defined by c(M) := ||M|a||, where M|a is the restriction of M to 
the subspace of vectors f * fulfilling YJi=i = 0. 

Clearly < c(M) < 1, see ^. Now, I may consider some special situations 
delivering more information about evident solutions constructing. 

Example 2. The most simple one is M{t) = M = const.. Then equation (|^ 
is the following: 

M = u{M -1)M. (13) 

Then for M = M"^ the solution exp(i^(M — l)t) has the limit for t ^ oo oi 
the form 

In the case c(M) < 1 the limit is O, where c(0) = 0. 
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Example 3. Let M{t) = a{t) ■ 1 + (1 — a{t)) ■ 0, where a, j3 are an integrable 
on finite intervals functions of values in the interval [0, 1]. Then equation 
is reduced to 



/3(T) = e-^^a(T) + i^e"''^ a(T - t)(3{t)e'''dt, (15) 

Jo 

where (3 is uniquely defined through M(t) = P(t) ■ 1 + (1 — P(t)) ■ O. In the 
similar way for each Mi, M2 G 2S such that {aMi + (1 - a)M2, a G [0, 1]} 
is closed for matrices multiplication linear integral equations arise of the 
form (3 = flaf^, where Qa transforms any integrable input function /3 of 
values in [0, 1] into an output QaP(t) G [0, 1]. The solutions of the equations 
(generically) stabilize in infinity to a number from [0, 1] as one may conclude 
from the asymptotic theorem in section 

In the case of locally constant functions a{t) except finite number of dis- 
continuities in bounded intervals equation ( [T3| ) may be viewed as a sequences 
of differential equations defining and solving step by step. I consider a simple 
0, 1 input function. 

Example 4- Let a{t) = a2k = 1 for t G [2kT, {2k + l)r) and a(t) = a2k+i = 



for t G [{2k + l)r, {2k + 2)r), where A; G N and r > 0. Then (|T5D is 
transformable to 

/3{T) = uj2(3{T- kT)e-''^^{ak - a^-i) (16) 

k=l 

with T G [zr, {i + l)r) or adopting periodicity 

/3(T + 2r) = e-2'^"/3(T) - vl3{T + r)e-"" + z/^(T)e-2^" (17) 

The initial conditions for each intervals are /5(0) = a(0) and the discontinuity 
in t = fcr are given by (3{kT^) — l3{kT~) = (— l)'^ exp(— z/A;r) for k > 1. In the 
way solving the integral equation in an interval one needs to possess already 
solutions of last two. First two intervals need to be solved independently. 
Here (3{T) = 1 for T G (0, r) and /3{T) = 1 + {ur - uT - 1) exp(-z/r) in 
[r,2r). 

For effective finding of solutions I return to an analytical input. 
Example 5. Let 

1 cos(t) 1 e** e~** , , 

a{t) = - + ^ = - + — + 18) 

^^22 244 ^ ' 
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and u = 1. 1 assume that the solution of (^) has a form 

e'p{t) = a{t) + b{t)e'' + 6(t)e-'* (19) 

where a, b smooth functions, a{t) G M, b{t) G C. Then the linear, constant 
coefficient equations follow 

a^^) - a + d - = 0, (20) 

b^^^ + {3i - l)b - 2{i + l)b - ^b = (21) 

with the initial conditions: 

a(0) = d(0) = d(0) = ^, (22) 

m = I, (23) 

m = (24) 

6(0) = 1(.-1). (25) 

Similarly, the equivalent, linear, constant coefficients differential system of 

equations of order 2A^+1 may be constructed for M{t) = Y.n=-N exp(mt), 
where An are constant matrices. 



IV Asymptotic theorem 

1 define a restriction on 2S'*+ maps 

Definition 2. Let M : M+ ^ M(n x n, M) be a measurable map into doubly 
stochastic matrices. It is said to be generic iff 3(5 < 1) fi^Is) > 0, where 
Is '■= {t G M+;c(M(t)) < 6} and fi the Lebesgue measure of the real, non- 
negative numbers space R+. 

The above genericity condition naturally appears for M{t) generated by 
physical quantum systems defined in section for details see appendix 0. 
Now, the main theorem may be proved. 

Theorem 1. Let M{t) be a generic map. Then limj_»oo M(t) = O. 
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Proof. One has the following estimation 



c{M{t)) < E / dtn-l ... dtiC{M{T - tn)) O... 

„=o-^o Jo Jo 

oc{M{t2 - ti)) o c(M(ti))z/"e-''^. (26) 



so the problem is reduced to considering equation (|T^) with a differing from 1 
on a positive measure set. Changing variables in the integral of the equation 
and making transformation T jv lead to elimination oi v (z/ = 1). Let 

/•oo 

5= / a{t)e-'dt. (27) 



From the assumption about a is that 6 < 1. Now, let (3{T) < bo for T G 
[to, oo). From (^) I have 



p{T) < e-^a{T) + e"'^ / " a{T - t)[3{t)e'dt + h^S. (28) 

Jo 

I define e = (1 — 5)/2. Let tg, > to, be such that two first terms of (^8|) are 
less or equal to boe for t > tg. Then /3(T) < 6o(^ + e) for T > tg. Continuing 
the procedure one reaches lim^^oo /5(7') =0. □ 

Nevertheless, even for c(M(t)) = 1 an asymptotic stabilization can exist 
as the example is showing. 

Example 6. Let M{t) = P, where P is a cyclic group generator such that 
^ 1 for 2 = 1, . . . , A; — 1 and P'' = 1. I also assume that c(X;(Li «i-P') = 1 
for ctj > and Y.i=i = 1. It may be easy realized by the matrix 



Pa 
P. 



(29) 



P„ is a permutation. Then the solution (^) has the form 

M{t) = exp{-ut)[P''Mut) + Ph{vt) + ... + P'~'Mut)]P, (30) 

where the analytical functions, fi{z), are formally defined by exp(aj2;) = 
(«,) Vi(z) + . . . + («0'"Vfc(^) with af = L Then lim^^^o M{t) = i Eti 

Therefore, a generalization may be proposed. 
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Theorem 2. Let M G 2E'^+ be a right-side continuous map and M(0) = 1. 
If M{t) fulfills the integral equation then 



lim M{t) 











o\ 





62 



















V 








id J 



(31) 



Proof. I put Sm '■= {M(t); t G Let Sm be the closure of Sm with respect 
to matrices muhiphcation and 3 := {Pa', ex G /} be a minimal permutations' 
subset spanning Sm- Then = 0j © Vid, where 7^ {0} are minimal 
invariant subspaces of P o Aj, P~^ G 3, spanned by the canonical basis 
vectors and such that dimVj > 1, where denotes all doubly stochastic 
matrices obtained from 3 as baricentric points, see proposition ^. On Vid, 
Vid -i- Vi, elements of acts as identity. 1 G A^, so P = 1 may be chosen. 
I define a compression in Vi by Cj(M) := c(M|vJ for M G Aj. If V(M G 
A^) Cj(M) = 1 then P' o Mlv^ is decomposable for P' being an admitted 
permutation of the basis vectors from Vi. Therefore, for all interior M G 
the compression Cj(M) < 1, also the next consequence is that 9j G A^, 
where I have denoted Aj := A:}\v^. Now, I restrict considerations to one 
subspace Vi and for simplicity omit the index. Let & C Sm be a minimal set 
that does not belong to any yl^, Z 3. & is finite and c{J2k (^kSk) < 1 for 
«fc > 0, J2k(^k = 1 and {Sfc} = &. Each has a form Mit^,.) o • • • oM(ti^). 
Let := max{Nk}. At the moment I return to equation (^ with u = 1. 
Equivalently, it may be written in the form: 



+e" 



-T 



M(T) = e 

T ftM l't2 



JO 



I will check the limit of W{T) : 
and 



MiT) + e"^ M(T - ti)M{ti)dti + ... (32) 

JO 

M(t2 - ti)M(ti)e*Mti . . . citTV- 
TniM{T + TO, where 7; > 0, Ei7i = 1 

(33) 



nk=l 



k are indices of 5^. 

The final step is like as in theorem |1|. Having an upper bound 60 of 
c{W{t)) in [to, 00) one improves it, here using (^2]), to (1 + 6)bo/2 on [tg, 00) 
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in the next step, where t'^ > to + 1 may be found. Now, I will show it. The 
matrix W{T) has the form: 

W{T) = J2 li[e~^^^^^^M{T + Ti) + f^^' M{T + Ti- t,)M{h)dh + ■■) 

+ E 7.^-^^+^') r^"' r ... r + - o . . . 

^ JO JO JO 

oM(t2 - ti)M(ti)e*Mti ...dtN + 
E 7.e-(^+^') Z^'"'^' r M{T + T,-t^)o... (34) 

^ JO JO Jto 

oM{t2 - ti)M{ti)e^'dti . . . dtN. 

One may verify that 
c(E li f M{t^^ - ejv) o . . . o M(ti, + £2 - ei)M(T + ei)e-^'+^^d^e) 

i JA 

<E7ze-'^'/x(A)supc(M(T + ei)), (35) 

where A is a measurable set of M^, /i(A) > and Si = Mit^i) o . . . o M(tiJ 
with additional M(0) = 1 if it is necessary for A^"; = A^. In the way I obtain 
an estimation for the last term of (|3^ ) by 6bo, 6 < 1. The first terms of (p^ 
have vanishing compression for T —>■ oo. The proof is completed. □ 

In other words the limit doubly stochastic matrix projects initial state p onto 
the maximal entropy state available for p. Theorem ^ covers all analytical 
situations arising from quantum physics. Nevertheless, for keeping the result 
( [HT| ) it is also enough to assume that M[t) differs from maps in theorem ^ 
on a set of measure zero. 

V Summary and interpretation 

I have modified the ordinary quantum evolution of a system-bath through 
considering an associated stochastic process describing probability transfor- 
mation, not Hilbert space vectors. The arising formalism offers an integral 
equation which transgresses the former methods of stochastic modification 
of differential equation based on adding a linear stochastic term. The closer 
analysis of the integral equation direct solving shows the difference with the 
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ordinary differential one, its non-Markovian character. The equation appears 
in a natural way from physical considerations and reproduces the time arrow 
- the missing feature of unitary evolution. The relation to the physical origin 
is stressed by a simple observation. Let M{t) be of the period 27i. If I change 
the period by a general r > then M'{t) corresponding to M'[t) = M{—t) 
is equal to M[—t), but, now, the integral equation is governed by the new 
u' such that 

Tu' = 2txv. (36) 

This is a theoretical suggestion for a universal relation between the period 
of quantum wave and the associated coefficient of stochastic reduction. The 
interpretation may be found even for cases excluded from the asymptotic the- 
orem Namely, if M{t) contains an identity sector arising from eigenvectors 
of hamiltonian, then in the sector the decay effect does not appear at all. In 
the way z/ > influences only states which during evolution cease to be pure. 
Appearance of invariant sectors is related to observables commuting with the 
hamiltonian, which are not also disturbed by the stochastic modification. 



A Genericity conditions 

The aim of the appendix is to study the maps M{t) from ]R_|. into 2E arising 
from quantum evolution. Firstly, I start from a characterization of doubly 
stochastic matrices with unit compression. 

Proposition 3. Let M he a doubly stochastic matrix. Then c(M) = 1 iff 
P o M is a decomposable matrix for P being a permutation. 

Proof. The implication (<^=) is obvious. If M" = © V2 is a decomposition 
then it is enough to put p = pi® qp\, g 7^ 1 to have p ^ p^ and ||M(p)|| = 
||PoM(p)|| = IIpII, where {p^^Y = l/dim\/ and (p^)^ = 1/diml^j, i = 1,2. 
To show (^) I take M = J^iCaPi, where = ca ^ and Pi are 

permutations enumerated by i. I may assume that ai ^ 0. Then Pf ^ oM = 
ail + Pi^ o CiiPi. Still c{P{^ o M) = 1, so p, p G A, exists such that 
P{^ o P^p = p and ||M(p)|| = Then the canonical basis vectors for 

which corresponding components p^ are equal to themselves constitute the 
invariant subspaces and a decomposition is done. □ 

From the above proposition one states that c(M) < 1 is open and dense in 
2S. 
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Doubly stochastic matrices can be also built via Kraus representation^. 
Then one begins from defined by equation (|5|). Generic maps 

M(t) defined by appear in the following way. 

Proposition 4. // at least one -B^^ G {Bab} is not proportional to the iden- 
tity then a dense and open set of orthonormal basises of C" exists such that 
M{A{t)) is generic. 

Proof. Now, M{t) = M{A(t)) defined by an orthonormal basis {\i >} ap- 
pears to be analytical, so of the form 

oo 

M{t) = E yM,, (37) 

fc=0 

where Mk G M{n x n, M). First terms of the expansion are Mq = 1, Mi = 
and 

(M2),7 = ^(E I < 3\Bab\l > I'- < j\BacBjl > 6,i). (38) 

a,b,c 

Basises {\i >} constituting the open and dense set of U{n) are such that 
j 7^ / ^ < j\B'^fj\l >7^ 0, where is indicated in the assumption. One 
finds that 

\\Mit)p\\ = 1 + Ep.(M2).,p//2 + ait'), (39) 
where = 1. Further 

j:MM2hP, = E E I < ^\Bab\j > r(m - < 0. (40) 

Therefore, for a generic basis {\i >} from ( pUD one obtains c(M(t)) < 1 for a 
t > 0. The genericity of M{t) is proved. □ 
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